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Competitive random sequential adsorption on a line from a binary mix of incident particles is 
studied using both an analytic recursive approach and Monte Carlo simulations. We find a strong 
correlation between the small and the large particle distributions so that while both partial con- 
tributions to the fill factor fiuctuate widely, the variance of the total fill factor remains relatively 
small. The variances of partial contributions themselves are quite different between the smaller and 
the larger particles, with the larger particle distribution being more correlated. The disparity in 
fiuctuations of partial fill factors increases with the particle size ratio. The additional variance in the 
partial contribution of smaller particle originates from the fluctuations in the size of gaps between 
larger particles. We discuss the implications of our results to semiconductor high-energy gamma 
detectors where the detector energy resolution is controlled by correlations in the cascade energy 
branching process. 
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I. INTRODUCTION 

One dimensional irreversible random sequential ad- 
sorption (RSA) has been of interest for several decades. 
Its numerous extensions include RSA with particles ex- 



panding in the adsorption_process 



tide adsorption 



two-size par- 



a III, and also RSA with an ar- 
bitrary particle-size distribution function j^. The inter- 
est is due to the relevance of this process to a number 
of physical phenomena in different fields of application, 
such as information processing 



impact ionization 



10 



under external stress 



9[, particle branching in 
and crack formations in crystals 
ll| . The simplest example of RSA 



is the so-called car parking problem (CPP). In the con- 
text of CPP, one studies the average number of particles 
( "cars" ) adsorbed on a long line and the variance of this 
number. Equivalently, one is concerned with the distri- 
bution function for the size of gaps between the parked 
cars (see Refs. [l^Q for the review). 

The problem of competitive RSA from a binary mix- 
ture is of special interest because of the non-trivial corre- 



Isubashiev® ece . sunysb . edu | 



lations in both the particle and gap-size distributions, de- 
veloped during the deposition. These correlations mani- 
fest themselves in the final irreversible state correspond- 
ing to the so-called jamming limif — when every gap 
capable of adsorbing a particle has done so. Numerous 
studies, reported in the literature for the binary-mixture 
RSA in the jamming limit, addressed the problem of cor- 
relations only indirectly, through its manifestation in the 
fill factor or the gap distribution. Available results in- 
clude binary mixtures with point-like particles [J, |5[ and 
those with a relatively small particle size ratio, b/a < 2 
J|. Also available are Monte-Carlo studies of the fill 
factor and the gap-size distribution for a binary-mixture 
deposition with equal abundance of both particles 

The present study is concerned with the correlation 
between the fluctuations in the number of adsorbed par- 
ticles of each kind from a two-size binary mixture, as 
well as with their partial contributions to the fill factor. 
We present both analytical results and those obtained 
by Monte-Carlo simulations for a wide range of binary- 
mixture compositions and size ratios. 

We are interested in the RSA problem primarily be- 
cause of its relevance to the propagation of high-energy 7- 



particles through a semiconductor crystal — with parti- 
cle energy branching (PEB) due to cascade multiplication 
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171 . The 



of secondary electrons and holes 
correlation of energy distribution between secondaries is 
quite similar to that of the gap distribution in the RSA 
process [l^. In both cases, the ratio of the variance of 
the final number of particles to the average particle num- 
ber in the final (jamming) state can be much less than 
unity, which is favorable for the detector energy resolu- 
tion. This ratio (which would be unity if the particle 
number obeyed a Poisson distribution) is called the Fano 



factor, $ 
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The reported attempts to evaluate <& employed over- 
simplified models of the semiconductor band structure. 
In such models, all crystal properties are characterized 
by three parameters, namely, the band gap, the phonon 
frequency, and the ratio of the rate of phonon emission to 
that of impact ionization. The price of this oversimpli- 
fication had been that correspondence with experiment 
could be achieved only by assuming unphysically large 
rates of phonon losses (about 0.5 eV per created e-h pair). 
This docs not corroborate with the known values for the 
ratio of the impact ionization and the phonon emission 
probabilities for high-energy electrons in semiconductors. 
The model furthermore obscures the role of features in 
the band structure and the ionization process that are 
specific to a particular semiconductor. 

In our earlier work ^ , we used an extended RSA model 
of particles that expand or shrink upon adsorption. The 
shrinking model is relevant to the PEB problem in that it 
helps to elucidate such factors as the non-constant den- 
sity of states in the semiconductor band and the fact that 
due to momentum conservation the ionization threshold 
is larger than the actual (bandgap) energy that is lost in 



impact ionization. 

The recursive technique employed in Ref. [3| allowed 
us to assess the accuracy of approximate approaches to 
the yield and variance calculations (such as, e.g., the 
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m 



average-loss approach of Refs. 

In the present work, the RSA model is extended in a 
different direction — competitive deposition of different- 
size particles from a binary mixture — that is suitable 
to simulate the role of multiple channels of pair produc- 
tion, owing to the multi-valley nature of semiconductor 
bands. We arrive at a number of qualitative conclusions 
that should be taken into account in both the interpre- 
tation of experimental data and the choice of the crystal 
composition and device structure in gamma detectors op- 
timized for energy resolution. 

The paper is organized as follows. Section II presents 
the basic equations of the recursive approach and the an- 
alytical results for the fill factor and its variance for the 
larger particles. In Sect. Ill, we analyze the results that 
demonstrate high correlation in the particle distribution. 
Based on the gained understanding, we formulate in Sec. 
IV the implications of our results for the Fano factor of 
semiconductor 7 detectors. Our conclusions are summa- 
rized in Sect. V. Certain analytical results are derived in 
the Appendix. 



II. PARTIAL CONTRIBUTIONS TO THE FILL 
FACTOR AND ITS VARIANCE FOR TWO-SIZE 
RSA PROBLEM 



We consider the problem of competitive deposition 
from a binary mixture of particles with sizes a and b, 
whose relative contributions to the total flux on the ad- 
sorbing line are q and p = I — q, respectively. We shall 
use a recursive approach to first study the mean number 



of particles na{x) and ni,{x), adsorbed on a line of length 
X (in the jamming limit), and then the corresponding 
variances. 

Consider a large enough empty length x > a,b. We 
assume that the adsorption is sequential, i.e. only one 
particle is adsorbed at a time. The first adsorbed par- 
ticle will be of size a with the probability of landing at 
any point q{x — a)/{x — I) or of size b with the landing 
probability p{x — b)/{x — I). Here I ~ qa + pb is the 
"average" particle size in the binary flux. After the first 
particle is adsorbed, it fills a certain interval [y, y + a\ (or 
[y,y + b]), and leaves two independent segments, whose 
combined size is either x ~ a or x — b. The average num- 
bers of a-particles Uaiy) and na{x — y — a) (or na{y) and 
na{x~y — b)) will be subsequently adsorbed in these gaps. 
Thus, the recursion relation is of the form 

na(x) = —[l + na(y)+na(x-a-y)\ 



X ~ I 

p{x — b) 
X — I 



na{y) + na{x - b - y)] 



where the first and the second terms (upper and lower 
lines) correspond to the cases of the first landed particle 
being a particle of sort a or 6, respectively. These cases 
must be averaged over all possible landing coordinates y 
of the first particle in a different way, viz. for a first, 

1 



< na{y) >a= 
whereas for b first, 

< na{y) >b= 



na{y)dy , 



1 



- b 



x — b 



na{y)dy . 



Performing the average and using the symmetry between 
left and right segments we obtain, finally: 
q{x — a) 2q 



na{x) 



X — I X — I Jo 

c-b 

na{y)dy. 



na{y)dy 



2p 



X - I Jo 



(1) 



A similar equation holds for the particles of size b: 
nb{x) = =- H = / nb{y)dy 



X — I X — I Jo 

/*x — h 



2p 



nb{y)dy. 



(2) 



x-l Jo 

With the help of Eqs. (|1I2[) one can readily derive an 
equation for the average total covered length f{x), de- 
fined as f{x) ~ ana(x) -f brn,(x), giving 

xl^qa^-pb^ 2q [^-^ 
fix) = = + = / .f{y)dy 

X ~ i X — I Jo 



2p 



I{y)dy- 



(3) 



Equation ^ agrees with that of Ref. [s] for the total cov- 
ered length in RSA from a multi-size mixture. However, 
the advantage of Eqs. (|ll2p is that they permit studying 
the partial contributions to the coverage by each of the 
two sorts of particles separately. 

Note that the symmetry between the a- and the b- par- 
ticles is broken by the initial conditions. To be specific, 
let b > a. Then, for 6-particles the boundary condition 
at small x is simply 



nb{x) = 0, < X <b 

whereas for a-particles we have 

0, < a; < a 



(4) 



'^"('^^ ^ ^ l[ a<x< min(2a,6) 



For b > 2a, Eq. ([5]) should be supplemented with 



na{x) = 1 + 



X — a 



na{y)dy 



(6) 







Eq. ^ accounts for the deposition of smaller particles in 
small gaps where the larger particle does not fit. Clearly, 
this process is not influenced by the 6-particles and does 
not involve particle competition. 

More refined arguments are needed to derive the sec- 
ond moment of the distribution, i.e. the expected value 



of the square of the number of particles of a given sort, 
Ua{x) = En^(a;). It may not be a priori evident that one 
can write independent expressions for particles of both 
sorts, because parameters a and b not only describe the 
particle size but also designate the sort of a particle. In- 
deed, we can even have a = b and distinguish the parti- 
cles by some other parameter, like "color" . Our approach 
should remain valid in this case too. To be rigorous, we 
therefore introduce an artificial parameter, the "mass" of 
a particle, and nib, whose value may depend on the 
particle shape and is simply proportional to the particle 
length only for a fixed transverse particle size. Hence one 
can regard and rrif, as independent parameters. 

Consider a total mass M{x) ~ mana{x) + m[,nb{x) 
of the particles adsorbed in a line segment x. We first 
evaluate recursively the mean square of the total mass 
< M'^{x) >= {[mana{x)+rrn,nb{x)Y) , and then calculate 
the second partial derivatives with respect to nia and mt,. 
Using the landing probabilities of particles to perform the 
averaging, we obtain 



Ua{x) = (x-l) 



^ q{x - a) + 2q / Ua{y)dy 
Jo 

j*x~h px—a 

+2p / Uaiy)dy + 4:q / na{y)dy 
Jo Jo 

+2q na{y)na{x - y ~ a)dy 

Jo 

px — b 

+ 2p na{y)na{x - y - b)dy 
Jo 



(7) 



Similarly, equation for Ub(x) reads 



Uh{x) = {x ~ ly 



p{x -b) + 2q / Ub{y)dy 



nx—b nx—b 

+2p / Ub{y)dy + Ap I nb{y)dy 
Jo Jo 

px—a 

+2q / nt{y)nb{x - y - a)dy 
Jo 

px — b 

+2p / nb{y)nb{x - y - b)dy 
Jo 



(8) 



We could have derived Eqs. (|ll2p in a similar way, 
by first evaluating the total average mass M{x) = 
"nianaix) + mbnb{x) recursively, and then calculating the 
derivatives. For a more general case, when the total mass 
is a linear functional M{x) = J mini{x)dl on the mass 
distribution mi , one would have to use variational deriva- 
tives dM{x)/6mi. For the case of binary mixtures we 
consider, partial derivatives are sufficient. 

Similarly, we derive an equation for the correlation 
fimction Uc{x) ~ {na{x)nb{x)) by calculating a mixed 
derivative of < > with respect to and nib- 

For particles uniform in the transverse direction with unit 
mass density, both the mass and the length of particles 
are identical, which gives a way to check the equations. 
An appropriate linear combination of equations for Ua, 
Mf,, and Uc then gives an equation for the variance of the 
total filled length or, equivalently, for the variance of the 
wasted length, w(x) — x — fix). The resulting equation 
can also be obtained directly, by applying recursion ar- 
guments to the waste. The identical results obtained can 
be viewed as an additional proof of Eqs. (|7I8|) . 

Note the asymmetry in the 4-th terms of Eqs. (|7l8p 
that are proportional, respectively, to 4g and Ap. These 
terms ensure the correct (linear) asymptotic behavior of 
the variance at large x. 

An important feature of Eqs. (|1I2I7I8P is that in spite 
of the competitive character of the deposition of particles 
of different sorts, the equations for ria, Ub and the higher 
moments are independent. This is rooted in the fact that 
a single deposition step on an empty length x does not 
depend on the already adsorbed particle distribution. 

Due to the self-averaging nature of the filling length 
(and waste length) in the limit a; ^ oo the averaged 
(hence approximate) recursion equations yield exact re- 
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FIG. 1: (color online) Average number of adsorbed particles 
Ha (solid lines) and ni, (dashed lines) as functions of the length 
X (measured in units of a) of the adsorption interval, assumed 
initially empty. The results are obtained by iterating Eqs. 
(|1I2|) with the assumed ratio of the particle size h/a = 2.4 
and the varying fraction q of a-particles in the flux. 

suits. The recursive technique is in this sense equiva- 
lent to the alternative "kinetic" approach to RSA that 
is sometimes regarded as a higher-level theory. In the 
kinetic approach one considers the rate equation that 
describes the sequential deposition of particles with the 
particle distribution on a line characterized by a time- 
dependent function G{x^ t) representing the average den- 
sity of gaps whose size is between x and x -\- dx 
It has been ascertained for a number of problems that 
both approaches give the same result for the coverage. 
Still, each has its own benefits. The kinetic approach 
allows studying the temporal variation of a state with 
specified particle distribution. The recursive approach, 
while simulating a simplified version of the kinetics, al- 
lows to study more complex effects, such as variance of 
the adsorbed particles of different size. 

Evaluation of na{x) and nb(x) is readily done by re- 
peated iterations of Eqs. (|1I2|) . going from the small to 
progressively larger lengths x. Results of the numerical 
recursion are shown in Fig. 1 for a particle size ratio 
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b/a = 2.4 and varying q. 

The noteworthy features of the functions na{x) and 
n{,{x) are (i) the step-like features a,t x = x ~ b 
(which are replicated with ever smaller amplitudes at 
X = na + mb, where n and m are integers), (ii) the dip 
in the number of small particles na{x) at x = b, which 
increases with p, and (iii) the reduction of ni, with in- 
creasing q. We also note that for all q the behavior of 
both na{x) and nh{x) becomes very close to linear al- 
ready at a; w 7. 

The asymptotic behavior of na{x) and ni,{x) at large 
X can be obtained by multiplying Eqs. pi2p by a; — I and 
taking the derivative with respect to x. The resulting 
differential equations are satisfied by linear functions of 
the form 

Ua = aaix + 1) - q, Ub ^ ai,{x + 1) - p, (9) 

where and are arbitrary constants. When correctly 
chosen (by matching to the recursive solution) these con- 
stants become the partial filling factors. After the match- 
ing is done, the total filled length in the asymptotic limit 
is given by f{x) = 6x + [6 ~ 1)1, where 9 = aua + bub is 
the specific coverage. It is worthwhile to stress that the 
value of the asymptotic solutions ^ consists precisely in 
that they are asymptotically exact. Hence they provide 
a sanity check on any solution we could have obtained by 
a numerical recursion up to moderate values of x. 
Similarly, Eqs. ([THll) yield the variances at large x, 

Ua-nl = ^aix + 1) - qp[l + (6 - a)aa]^, (10a) 
Ub-nl = i^Lb{x + 1) - qp[l - (b - a)ab]'^. (10b) 

Again, these solutions are asymptotically exact; they sat- 
isfy Eqs. (|7l8p with arbitrary values of /ia and /it, pro- 
vided of course that na{x) and nb{x) are in the correct 



asymptotic form ([9]) with properly chosen [i.e., satisfying 
Eqs. (jll2p ] coefficients aa and at,. In principle, we could 
now follow a procedure similar to above, viz. determine 
/ia and /ib by matching Eqs. (jlOp against a numerical 
recursive solution at some moderate value of x. How- 
ever, it would be rather difficult to control the numerical 
accuracy in this procedure, because of the difference of 
nonlinear functions that enter Eqs. (jlOp . even though 
that difference itself behaves linearly with x at large x. 

Fortunately, our model admits of an exact solution 
based on the use of Laplace transformation (details can 
be found in [3[ and references therein). Below we present 
an exact evaluation of variance for particles of larger size, 
while details of similar though lengthier calculations for 
smaller particles are presented in Appendix. 

Firstly, we need exact solutions of Eqs. (jl|2p . To ob- 
tain these, we substitute a; — > 2; + in Eq. ([2]) and mul- 
tiply it by X — I. Taking the Laplace transformation of 
the resulting equation and using the boundary condition 
(III), we obtain 



h)/{h — I) at small x — b. Hence we have 



-— + b-l 



(11) 



Here Nb(s) is the Laplace transform of nif{x)^ 

m{s) = / e-''=nb{x)dx, (12) 



Rearranging the terms and multiplying by e we put 
Eq. pT|) into the form 



I + - (qe 



N,{s)^^^e~''\ (13) 



For p 00, the solution of Eq. (|13p is, asymptotically, 



Nb{s)\s 



s(&-0 

as follows from the known variation of ni,(x) « p(^x 



(14) 



Nb{s) 

where 
l3{s) = exp 



pcxp(— ? s) 
s2/3(s) 



I3{t)dt, 



1 — qcxjp{~at) — pcxp(— 6t) 



t 



(15) 



dt . 
(16) 



To find the asymptotic behavior of ni,{x) at large x, it 
is convenient to use Karamata's Tauberian theorem for 
the asymptotic growth rate of steadily growing functions 
(see e.g. 



20[, p. 37). According to the theorem, the 
asymptotics of n(,(a;) [or na{x) or their variances] can be 
readily obtained (by taking the inverse Laplace transfor- 
mation) from the Laurent power series expansion of the 
Laplace transforms of these functions at small s (see Q 
for the mathematical details of this analysis). 

Function Ni,{s) is analytic at all s and at s = it 
has a second-order pole with the following asymptotic 



Nb[s) = H h 0(s) 



where 



"6,0 = P e 
Iq 



— q{b—a)s 



P{s)ds. 



(17) 



(18) 



To calculate nb{x) at large x, we take the inverse Laplace 
transformation of P?|) . This gives 



nb{x) = abfl{x + 1) -p, 



(19) 



with an exponentially small error term, in line with the 
asymptotics given by Eq. ([9|. 

In the limit p = 1, equation ((T8)) duly gives the so- 
called jamming filling factor R for the standard RSA, 
abfi{l = 1) = ^ = 0.74759- •• (also called the Renyi 



constant 



2l|). In the limit a ^ 0, Eq. p8)) recovers the 
results of Refs. P, 0] for the coverage of a line from a 
binary mixture of finite size particles and point defects. 



Moreover, Eq. (|18p gives the large particle contribution 
to the total coverage, obtained in 6, ^ for the range 
a < b < 2a. Here we see that this result remains valid 
for arbitrary a < b. 

Next, we perform similar manipulations with Eq. ([5]) 
and obtain an equation for the Laplace transform of the 
variance Ub{s) ~ L[uii(x)], viz. 



where 



I + - (ge""" +pe 
s ^ 



exp(-bs) 
Ub[s) = Rb[s), 

(20) 



Rb{s) =p + ipsNbis) + 2s^N^is) (geC'-'^)-^ +p) , (21) 

with iV6(s) defined by Eq. ([15]). The solution of Eq. 1^ 
can be written in a form similar to Eq. (|15p . namely 

Ub{s) = ^^f^ 1^ me-'^'-''^'Rb{t)dt. (22) 

The integrand in the right-hand side of Eq. ([22|l is pro- 
portional to l/t"^ causing the integral to diverge as 1/s 
for s — > 0. This is due to the square-law dependence of 
u{x) at large x . 

To separate the regular part needed for the estima- 
tion of variance, we note that at small t one has Rb{t) oc 
ctb.ot~^- Moreover, the series expansion shows that the 
difference P{t) exp[— ^(fe — a)t]Rb{t) — 2al ^t^^ is regular 
at t — > 0. Therefore, it is convenient to define an en- 
tire function Kb{t) ~ (3{t) cxp[—q{b — a)t]Rb{t) — 2al Qt~^. 
In terms of this function, the solution Ub{s) can be ex- 
pressed as follows: 
cxp(— si) 



Ub{s) 



where 



s2/3(s) 



kb,( 



+ h.o 



Kb(t)dt. 



Kb{t)dt 



(23) 



(24) 



To apply Karamata's Tauberian theorem, we note that 
the asymptotic expansion of Ub{s) near its third-order 



pole is of the form 

2ab,o f^b,o + 2a^ gl 



Ubis) 



kb.ol - Kfc(O) - qp(b - a)'^al Q 



(25) 



Taking the inverse Laplace transformation, we find the 
asymptotic form of Ub{x): 

Ub{x) = al^^x^ + {kbfi + 2alQi)x 

+kbfil - Kb{Q) - qp{b - afal^, (26) 



with an exponentially small error term. Using Eq. ([T^ 
to subtract n1{x). we find an equation of the form (jlOb[) 
with pb = kbfl + 2pabfl. The specific variance of the 
adsorbed number of 6-particlcs is given by (at x — > oo) 



Hb = a,,,o(l + 2p) + 2 / {/3(s)siVfc(s)e'^ [2pe 



-bs 



+sNb{s) (ge-"^ + pe"^'^)] -'^\ds. (27) 



Integrating by parts the last term and rearranging the 
result, we finally obtain 

Mb = «b,o(l - 2p) + 4p/„°° H^LMe-f" (1 - qe-- 

-pe-»-) du + 2 ^e-'-K{u)du, (28) 

where 



-an _p^~bu^ _ {^a + l)u\ 



K{u) = qe"''" [2(1 - qe' 

+pe-^'' [2(1 - ge-"" - pe"'"') - (6 -t- l)u] (29) 

and 

ab{u) = abfl-p I e-''^'-''^ypiy)dy. (30) 

In the limit of small p ^ 0, the Fano factor <& = 
fJ'b/ctb.o ^ 1- In this limit, large particles are distributed 
on the line randomly, without correlations. In the oppo- 
site limit, p = 1, Eq. (p8)) reduces to the standard RSA 
result, first obtained for a lattice RSA model by Macken- 



221. The numerical value of the Mackenzie constant. 
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fio = 0.0381564- •■, corresponds to $ = 0.0510387 ■•■ , 
see Expression ([^5]) for the larger particles has the 
same structure as the corresponding formula in the stan- 
dard RSA model (fixed-size CPP). Due to the exponen- 
tial factors in the integrands of Eq. , the dependence 
of fib on a for a <C 6 is quite weak. The limiting value of 
the specific variance for a/6 gives the specific vari- 
ance of the fill factor for the case of finite-size particles 
(6 = 1) mixed with point-size particles, 

^^b,p = a6,p,o(l - 2p) + V fiiL^e-" (1 - e"«) du 
+2pr|dS{'/e--(2-2e-«-.) 
+e-(i+P)" [2p (1 - e"") - (1 + p)u] } du, (31) 

where ab^p^ is the fill factor for this case, 



ab,piu)^p e /3p{y)dy, ab,p,o ^ ab,p{u ^ 0) 



(32) 



and 



Pp{u) = exp 



(33) 



It is worth to note that Eqs. (|2l8p and their solutions can 
be readily generalized to the case when particles of the 
smaller size have an arbit rary distribution in the interval 



[oi , 02] so long as 02 < 6 



23| 



The above analytic results for the variance of larger 
particles are essentially exact, as will be confirmed in the 
next Section by Monte Carlo simulations. For the smaller 
particles, the calculations are messier and accurate ana- 
lytical results can be obtained only in a certain range of 
particle size ratios. Estimations of the variance for small- 
size particles are further discussed in the Appendix. 



III. DISCUSSION OF THE RESULTS, 
COMPARISON WITH MONTE CARLO 
MODELING 



Here we present the results of numerical calculations 
using both the analytical expressions obtained in the pre- 
ceding section and Monte Carlo simulations. For large- 
size particles the Monte Carlo results are very close to 
analytical expressions both for the fill factor and the vari- 
ance, so we shall not dwell on their comparison. For 
small-size particles, especially in the range 2 < 6/a < 8, 
analytical calculations are rather unwieldy, so Monte 
Carlo simulations become indispensable. Larger size ra- 
tios, b/a > 8, lend themselves to an approximate ana- 
lytical approach (see Appendix). In this case, we use 
the Monte Carlo to estimate its accuracy for the small 
particle contribution. 

Traditional studies of the generalized RSA via Monte 
Carlo simulations follow a temporal sequence of events. 
For the case of adsorption on a line of the length x from 
a binary mixture, one step of the sequence comprises: 

(i) selection of a particle from the mixture according 
to the deposition flux ratio (with the probability q of 
choosing the small-size (a) particle, and the probability 
p = 1 — q oi selecting a particle of larger size b) ; 

(ii) random choice of a deposition coordinate of particle 
center on the line x with formerly deposited particles; 

(iii) rejection of the particle if it overlaps by any part 
with formerly deposited particles or with the line bor- 
ders; otherwise, the particle deposition proceeds with the 
formation of two new disconnected adsorption lengths. 

This traditional approach has several drawbacks, that 
make the modeling very demanding, both in terms of the 
computer time and memory allocation. 

Firstly, both the flUcd length in the jamming limit and 



the specific fill factor (coverage) depend on the initial 
length. Due to the self-averaging property of the coverage 
it tends to a unique exact value in the limit a; ^ oo. To 
obtain the accuracy of about 0.1 %, the common strategy 
has been to use large initial length values (10^6 -10^5) and 
make additional averaging over a set of about N^^ =100- 
1000 different realizations. 

Secondly, as time evolves and the jamming limit is ap- 
proached, the probability of finding a free gap for parti- 
cle deposition becomes greatly reduced, so that the ad- 
sorption time tends to infinity. The process is termi- 
nated when variations of the adsorbed particle number 
are smaller than those required by the desired accuracy. 

The recursive analysis of the generalized RSA suggests 
a revision of the above scheme. Since the deposition is 
random and sequential, it does not depend on the tem- 
poral history of the process or the growing number of 
rejected particles and their coordinates. Therefore one 
step of the sequence can be chosen as follows: 

(i) selection of any free deposition length, li > a. It is 
convenient to choose for li the outermost free deposition 
length on the left-hand side. 

(ii) if li < b, then particle of size a is deposited, 
otherwise the deposited particle is chosen according to 
the landing probability, given by q{li — a)/{li — I) for 
a-particle and p{li — b)/{li — I) for ^-particle, where 
I = qa + ph. 

(iii) random choice of a deposition coordinate (taken 
as the coordinate of particle's left end) on the line h for a 
given particle size, i.e. within the interval Zi — a for a-size 
or within li — b for 6-size particle, with the formation of 
two new adsorption lengths from the initial length li. 

It is readily seen that although the sequence of deposi- 
tion events is different from the actual temporal sequence 
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of adsorption (the simulated deposition proceeds by se- 
quentially filling the left-hand lengths), the statistics of 
divisions is identical and therefore so is the final distri- 
bution of the gaps, as well as all statistical properties 
of the jamming state. Our sequential scheme excludes 
deposition of to-be-rejected particles and therefore is in- 
comparably faster. Besides, it terminates exactly when 
the jamming limit (with no gaps larger than unity) is 
achieved. Direct comparison with the traditional Monte 
Carlo results, e.g. 0,3,01 exhibits total agreement. The 
difference in the calculation time is especially evident for 
small (close to zero) q: in the time scale of "real" deposi- 
tion, the jamming limit will be strongly delayed because 
of the rarity of events with small particle chosen. In our 
modified approach, all gaps smaller than b are "rapidly" 
populated by small-size particles, however small be the 
value of q. 

The next step of the revision is to exploit the fact 
(proven analytically in the preceding section) that in the 
jamming limit, the linear dependence on the adsorption 
length of both the average filled length and its variance is 
exponentially accurate, starting from a reasonably short 
length, certainly not exceeding x « 106. Since this lin- 
ear dependence has only two parameters [actually only 
one, as the parameter ratio is exactly fixed by analyti- 
cal considerations, Eqs. (|9|10p ]. both the coverage and 
the variance can be determined with Monte Carlo simu- 
lations of short samples. 

To be sure, in order to achieve the same accuracy as 
that obtained for long samples, the results should be av- 
eraged over a sufficient number of realizations iVj^. This, 
however, takes little memory or time. Calculations show 
similar accuracy for different x and Nfj, so long as their 
product X X Nfi is fixed. The results presented below were 
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FIG. 2: (color online) Partial contributions of small and large 
particles to the total coverage depicted as functions of q, for 
different particle-size ratios 6/a in the flow. Open points cor- 
respond to the limit b/a 3> 1, as described by the analytical 
formulae (|32p for 6-particles and (|A13|I for a-particles. For the 
total coverage, the open circles to the wasted length product 
approximation, Eq. p4p . 



obtained using a sample of size x = 200a for 6/a < 10 
and X = 400a for h/a ~ 20,40, subsequently averaged 
over 10 000 realizations, which appeared to be sufficient 
to eliminate any spread of the results in the graphical pre- 
sentation (producing an accuracy of better than 0.1%). 

The use of small samples is very effective in reducing 
the calculation time (with an ordinary PC, high-accuracy 
results can be obtained in minutes, compared to days in 



the traditional scheme 



Figure 2 shows partial contributions to the coverage as 
functions of the fraction q of small particles in the binary 
mixture at different ratios of particle size. As q increases, 
the coverage with large particles is substituted by that 
with small particles, producing some decrease in the to- 
tal coverage. In the regions of corresponding parameters, 
our results reproduce those of reported analytical calcu- 
lations (i.e. for b/a < 2 [^,0] and for a ^ for the 
large particle contribution) and those obtained by the 
Monte Carlo simulations of Q, Q, 0] , demonstrating the 
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FIG. 3: (color online) Variance of the partial coverage by ad- 
sorbed fe-particles from the binary mixture for different values 
of the particle size ratio b/a in the flow. 



validity of our revised approach. 

It is evident from Fig. 2, that the total coverage in- 
creases at smaller q, as can be explained by sequential 
deposition of the two kinds of particles. In the regime of 
small q, large particles are adsorbed first and their de- 
position, unobstructed by small particles, is tight. Sub- 
sequently, the small particle fill the gaps between large 
particles and this clearly reduces the total wasted length. 

The effect of increasing the particle size ratio b/a is 
pronounced only for b/a < 10, then it rapidly saturates. 
Therefore, for large b/a, say b/a = 20 the coverage by 
large particles is very close to that obtained for a model 
mixture of point-like and finite-size particles [by formally 
letting a = in Eq. (fT8|) ]. Such a model, however, has lit- 
tle relevance to any practical situation, because it simply 
ignores the partial contribution of small particles to the 
total coverage. The latter can be described analytically 
in the limiting case b/a —f oo, Eq. (jA13[) . 

The partial contribution of small particles steadily 
grows with the increasing size ratio due to the expand- 
ing gaps between the large particles. In the limit q ^ 0, 
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the total coverage can be estimated by observing that 
the specific wasted length in this case is a simple prod- 
uct of the specific lengths wasted in initial deposition of 
large particles and subsequent deposition of small par- 
ticles, i.e. 1 - 61 = (1 - 6a){l - Ob). Since for g = 
the specific coverage 9b ^ R and since for large size ra- 
tios (when the gaps between large particles are large) the 
specific coverage Oa = i?, we have 6* = 1 — (1 — i?)^=0.936, 
in agreement with the results reported in the literature 
y, Q]. However, the sequential nature of the deposition 
suggests that the entire q dependence of the total can 
also be approximated by a product of the specific wasted 
lengths in the competitive deposition of large particles 
q{\ — R) and subsequent deposition of small particles in 
the remaining gaps, which gives 

61 = 1 - (1 - i?)(l -pR) = R[l+p{l~R)]. (34) 

This product-waste approximation is shown in Fig. 2 by 
the open circles. 

Next, we concentrate on the specific waste variance and 
the Fano factor. We shall discuss the b- and a-particles 
separately, since the effects are rather different in na- 
ture and also since they have been evaluated by differ- 
ent techniques. Results for large particles are obtained 
by numerical integration of Eq. (|30p and confirmed by 
Monte Carlo simulations. Results for a particles are ob- 
tained by Monte Carlo stimulations and arc accompanied 
by analytical expressions in the limit b/a ^ 1. 

Variance, /2, of the partial contribution of 6-particles 
to the total coverage is shown in Fig. 3 for different par- 
ticle size ratios. Unlike the particle number variance ^, 
the variance of coverage, p. = fib, depends only on the 
size ratio b/a and does not directly scale with b. It is 
therefore more indicative of the effect of decreasing size 
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FIG. 4: (color online) The Fano factor for the number of ad- 
sorbed 6-particles from the binary mixture as functions of q, 
for different values of the ratio b/a of particle sizes in the 
flow. Dotted line corresponds to a random particle packing 
on a lattice with suitable lattice constant {aka monomer ad- 
sorption) . 

of small particles on the fluctuations of the number of 
large particles. At g — ^ 0, when the adsorption of large 
particles is unconstrained by small particles, the vari- 
ance of large particles is minimal and corresponds to the 
highly correlated distribution |l^ in the standard CPP 
problem (one-size RSA). The variance rapidly increases 
with q as the small particle deposition destroys the CPP 
correlations. The maximum of this effect is shifted to 
larger q values for larger b/a. For q approaching unity, 
the variance decreases simply due to the decrease of the 
average number of adsorbed 6-particles. 

Correlation effects are more adequately characterized 
by the Fano factor <&;,, shown in Fig. 4. With the increas- 
ing number of competing small particles in the flux, the 
Fano factor grows from the smallest value $ = 0.051 • • • , 
corresponding to the one-size RSA problem, to unity in 
the limit q —> 1. Small coverage by the large particles in 
the latter limit means that they arc distributed randomly 
on the line, so that Poisson statistics recovers. The most 
noticeable effect is a rapid decrease of the Fano factor 
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Small particle contribution to the flux, q 

FIG. 5: (color online) The Fano factor for adsorbed a-particles 
from the binary mixture for different values of the particle size 
ratio b/a in the flow. Open points show the contribution of 
fluctuations of the gap sizes between large particles. 

with 1 — manifesting a strong enhancement of the cor- 
relation effects in the large particle distribution. These 
correlation effects become exhausted only near q < 0.1. 
The correlation effects increase with b/a but saturate at 
about b/a = 20. 

Figure 5 shows the Fano factor for a-particles com- 
petitively deposited along with large particles. The re- 
sults are strikingly different at &\\ q ^ 1 (when $a = 
as expected). While the distribution remains correlated 
((f>a < 1) for small ratios b/a < 5, at larger b/a one has 
<i>a > 1, almost for all q, which means that the number 
of small particles per unit length is strongly fluctuating. 
This is due to the widely fluctuating size of the gaps 
available for small particle deposition between large par- 
ticles. For large values of b/a and in the entire range 
of q, the Fano factor $a can be approximated in terms 
of the fluctuations of the coverage by the large particles, 
viz. (f>a ^ {b/a)^b,pR'^/da, where fib,p is given by Eq. 
([31]) and 9a by Eq. (jA13p . This approximation, which 
neglects fluctuations of the density of adsorbed a parti- 




FIG. 6: (color online) Variance of the partial number of ad- 
sorbed a- and 6-particles and of the total number of adsorbed 
particles for b/a = 1.2. Also shown is the fluctuation correla- 
tion function fcor- 

cles in the gaps, is shown in Fig. 5 by open points. This 
contribution is proportional to b/a and for b/a > 10 it is 
evidently dominant. 

For the particle energy branching process at small 
b/a < 2, both the variance of the partial numbers of small 
and large particles and the total number variance are of 
importance. We shall illustrate this point in the instance 
oi b/a ^ 1.2 shown in Fig. 6. We see that a,t q ^ 0.5 the 
fill factor fluctuations are larger for a particles and some- 
what smaller for b particles, but both are pretty large, 
compared to the variance of the total number of adsorbed 
particles. This is due to the strong anti-correlation in 
their distribution, as evidenced by the specific fluctua- 
tion correlation function, fcor = x^^ (SnaSni,) , also plot- 
ted in Fig. 6. We note that fcor < 0, which means that 
any excess in the number of a-particles is accompanied 
by a downward fluctuation in the number of adsorbed 
6-particles. Importantly, the variance and the Fano fac- 
tor for the total number of adsorbed particles does not 
exceed substantially its value for the single-size RSA. 

Note the asymmetry of the curves for a and b parti- 
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cles, e.g. the variance of large particles goes to zero as 
q 1 whereas that of small particles remains finite even 
as g ^ 0. This is a feature of our model that allows "infi- 
nite" amount of time for the deposition of small particles 
in the gaps left after the deposition of large particles is 
completed, but not vice versa. Therefore, the deposition 
of small particles remains finite even in the limit of q — > 
and the same is true for the a-particlc number variance. 

Another interesting feature of the a-particlc num- 
ber variance, already evident from Fig. 5, is its non- 
monotonic behavior as function of b/a at small q. This 
variation is displayed directly in Fig. 7 that shows the 
dependence of the Fano factor on b/a for q=0.Q5, 0.1 
and 0.2 — where its non-monotonic nature is most pro- 
nounced. The minimum of the Fano factor is achieved at 
6/a « 2. Note that the non-monotonic dependence of the 
Fano factor is accompanied by non-monotonic variations 
in the dispersion of the gaps between small particles. In 
Ref. [7[ it was found that for q =0.5 the dispersion is 
noticeably reduced at b/a « 1.55. These effects were 
interpreted as a manifestation of the so-called "snug fit" 
events, i.e. particle deposition in gaps that are just barely 
above the unit length a. In contrast, the Fano factor for 
6-particles and that for the total number of particles re- 
main monotonic everywhere. 



IV. SOME CONSEQUENCES FOR THE 
ENERGY BRANCHING IN HIGH-ENERGY 
PARTICLE DETECTORS 



The model of RSA from binary mixtures is relevant 
to an important practical problem of particle energy 
branching (FEB) where high-energy particle propagates 
in an absorbing medium and multiplies producing sec- 
ondary clcctron-holc (e-h) pairs. Multiplication proceeds 
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FIG. 7; (color online) The Fano factor for adsorbed a- and 
fe-particles as functions of the particle size ratio b/a. Also 
shown is the Fano factor for the total number of adsorbed 
particles. 



so long as the particle energy is above the impact ioniza- 
tion threshold [ist . The energy distribution of secondary 
particles is random to a good approximation. 

The affinity between the two problems was fully rec- 
ognized already in 1965 by van Roosbroek |l^ (see also 
2J|). The FEB process can be considered in terms of 
a CPP if one identifies the initial particle kinetic energy 
with an available parking length and the pair creation en- 
ergy with the car size. Similarly, the kinetic energies of 
secondary particles can be identified with the new gaps 
created after deposition of a particle. Full equivalence 
of FEB to CPF further requires that only one of the 
secondary particles takes on a significant energy, which 



25[. Otherwise, one has 



corresponds to binary cascades 
to consider a simultaneous random parking of two cars 
in one event. 

To estimate the particle initial energy in FEB, one 
measures the number N of created electron- hole pairs. 
Variance of this number, due to the random character of 
energy branching and also due to random energy losses 
in phonon emission, limits the accuracy of energy mea- 
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surcmcnts. Both the yield TV and the c-h pair variance 



var(iV) = {N — N)^ are proportional to the initial en- 
ergy. The ratio of the e-h pair variance to the yield, i.e. 
the Fano factor of the FEB process, is a parameter that 
quantifies the energy resolution of high-energy particle 
detectors. 

For semiconductor crystals, the FEB problem has ad- 
ditional complications due to the energy dependence of 
phonon losses and the energy dependence of the electron 
density of states and the impact ionization matrix ele- 
ment. Full quantitative analysis of the FEB is possible 
only with detailed numerical calculations, which goes far 
beyond the scope of the present article. 

A common feature of the energy branching process in 
semiconductors is the presence of several pair production 
channels, associated with the multi- valley energy band 
structure of the crystal. In Si, Ge and common A3B5 
semiconductors, the e-h pair creation produces electrons 
in one of the ellipsoids near the edge of the Brillouin zone, 
in 100 (X) or 111 (L) directions. Owing to the difference 
in the final densities of states and the matrix elements, 
the impact ionization processes associated with X and L 
valleys have different but competitive probabilities. Be- 
cause of its low density of states, the F valley is usually 
not competitive, even when it is the lowest valley. 

Ultimately, electrons will end up in the lowest energy 
valley but when the final electron valley is itself degener- 
ate, as in Ge or Si, the resulting electron states may not 
be fully equivalent, because of the different collection ki- 
netics owing to the crystal anisotropy. This effect may 
have important consequences for the observed variance. 
For example, in Si diode detectors electrons are created 
in 6 degenerate energy valleys that represent ellipsoids 
of revolution elongated along (100) and equivalent direc- 



tions in k-space. Suppose the diode structure is such that 
the current flows along the (100) direction, as it is usually 
the case. Electrons from the two valleys along the cur- 
rent have a large mass and low mobility. The measured 
current is hence dominated by electrons from the 4 val- 
leys elongated perpendicular to the current that have a 
low mass and high mobility along the current. Since the 
choice of equivalent valley in the FEB process is fully ran- 
dom, the number of high-mobility electrons will fluctuate 
more strongly than the total number of generated carri- 
ers. These fluctuations will dominate if the inter-valley 
transition rate is low compared to the inverse collection 
time. In the opposite limit of high inter-valley transition 
rates, this effect will average out as the collected current 
will fluctuate in time. The current fluctuation mecha- 
nism due to the carrier escape into heavy-mass valleys is 
a well-known source of noise in multi- valley semiconduc- 
tors [26|). More detailed account for these effects will be 



presented elsewhere 
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Here we shall discuss an opposite situation, that is 
common to direct-gap semiconductors, such as GaAs or 
InF. In these materials, the lowest (F) electron valley has 
a very low density of states, compared to that in the 
satellite (X and L) valleys. Therefore, the probability of 
electron generation in the F-valley can be neglected in 
first approximation, so that the branching competition 
occurs only between the satellite valleys of two different 
kinds. Both the density of states and the threshold en- 
ergy are different between X and L valleys and we can use 
the results of the present study to interpret and predict 
the consequences, at least qualitatively. 

The binary-mixture RSA model interprets the higher 
density of states as higher deposition rate and the higher 
threshold as larger particle size. To make our conclu- 
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FIG. 8: (color online) Partial fill factors and the total coverage 
for b/a — 1.4 as a function of q. Also shown is the total 
number Ntot of adsorbed particles 



sions more transparent, let us re-formulate the required 
results in terms of a random parking problem with cars 
of two sizes. We are now interested only in the numbers 
of parked cars and the fluctuations of these numbers. 

Several qualitative conclusions can be drawn from our 
results: 

(i) The total number of parked cars (in the jamming 
state) will decrease with increasing fraction of larger 
cars in the flow and with the growth of their size. For 
b/a = 1.4 the effect is illustrated in Fig. 8 (which can 
be viewed as an extension of Fig. 2). It follows from the 
fact that adsorption of a large car excludes larger length 
for subsequent parking events and thus causes a decrease 
of the total fill factor. Note that the decrease in the to- 
tal particle number is accompanied by an increase in the 
total filled length, as smaller number of cars cover larger 
area. 

The next two conclusions (ii) and (iii), illustrated in 
Fig. 6, are interconnected and will be discussed jointly. 

(ii) Variance of the total number of parked cars and the 
Fano factor will both grow with the increasing fraction of 



larger cars in the flow and with the growth of their size. 

(iii) Variance of the separate numbers of parked small 
and large cars and their Fano factors are considerably 
larger than that of the total number of cars. Therefore, 
if for some reason one type of cars is neglected or under- 
counted, the registered variance and the Fano factor can 
be substantially increased. 

These conclusions are connected with the nature of the 
car number fluctuations and the strong anti-correlation 
between the fluctuations in the number of small and large 
cars. Fluctuations in the number of parked cars of one 
kind are strongly enhanced by the presence of more or 
less randomly distributed cars of the second kind, espe- 
cially when cars of the second kind dominate. This leads 
to conclusion (iii). However, the two distributions are 
anti-correlated (higher number of parked small cars is 
accompanied by a smaller number of large cars and vice 
versa). The anti-correlation is particularly strong for a 
size ratio that is close to unity. 

One can imagine a case when the two kinds of cars dif- 
fer only in "color" . In this case, Eqs. (|18l28p yield aa.o ~ 
qR, abfl = pR, i-ia = Rqp + q^fJ-o and = Rqp+P^fJ-o, so 
that at large x we have < SriaSni, > / x = —{R — fJ'o)'}P- 
Then, the anti-correlation is almost complete: the fluc- 
tuations of the total number are much smaller than those 
of a given color, but still non-zero. Both the individual- 
color number fluctuations and the anti-correlation are 
largest at q « 0.5, cf Fig. 6. The anti-correlation 
decreases with increasing size ratio, as reflected in our 
conclusion (ii). 

To discuss the above conclusions in terms of the FEB 
problem, we note that estimation of the initial particle 
energy is equivalent in CPP to a measurement of the 
unknown length of a parking lot in terms of the total 



number of cars that were able to fit into it by random 
parking, assuming tliat the average fill factor for a given 
two-size car mixture is known from earlier measurements. 
The absolute accuracy of such a measurement depends on 
the variance of the fill factor, and the relative accuracy 
is determined by the Fano factor. As shown above for a 
mixture of cars, the larger disparity of car sizes leads to 
the higher fill-factor variance and therefore reduces the 
absolute accuracy. 

A particle detector measures the total number of sec- 
ondary particles of all sorts (but not their total creation 
energy, that would be equivalent to the filled length). In 
any channel, all secondaries that have sufficient energy 
for further branching will do so. Therefore, only those 
pair creation energy ratios that leave the channels com- 
petitive (i.e. h/a < 2) are relevant to the FEB problem 
— otherwise additional energy branching would be pos- 
sible. 

We conclude that the presence of competing channels 
with different energies [e.g. impact ionization with ex- 
citation in X and L valleys] will decrease the quantum 
yield (the number of secondaries per unit energy of the 
primary particle) and enlarge the Fano factor. The at- 
tendant loss in energy resolution is not that bad when 
the ionization energies associated with different valleys 
arc not too disparate. For example, in Ge besides the 
lowest eight L valleys [Eq = 0.66 eV) one has a non- 
competitive r valley {Ey ~ 0.8 eV) and six very compet- 
itive Si-like valleys {Eq = 0.85 eV). The downgrading of 
energy resolution should be more important for crystals 
with larger (w 2) threshold energy ratio. For example, in 
Si one has besides the 6 lowest valleys {Eg = 1.12 eV) in 
X direction, eight germanium-like L valleys with the gap 
El = 2.0 eV. Their effect on the Fano factor in silicon 
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may not be negligible. 

Finally, reformulating (iii), we stress that any signif- 
icant disparity in the collection efficiency between dif- 
ferent equivalent valleys will strongly enhance the Fano 
factor and downgrade energy resolution. This happens 
because any collection disparity breaks the symmetry 
between the equivalent valleys and destroys the anti- 
correlation, responsible for keeping the total Fano factor 
low even when the partial particle numbers associated 
with individual valleys exhibit fully random fluctuations. 
One possible origin for the asymmetry in the collection 
efficiency in semiconductors has been discussed above in 
the case of silicon diodes with the electric field in (100) 
direction. In germanium diodes all different valleys are 
equivalent relative to the (100) direction and the sym- 
metry is not broken. It would be broken, however, if one 
were to use Ge diodes oriented in (111) direction. This 
would lead to a situation similar to Si — with a possi- 
ble degradation in the Fano factor. These effects deserve 
additional study, both experimental and theoretical. 

V. CONCLUSIONS 

We have studied a generalized 1-dimensional competi- 
tive random sequential adsorption problem from a binary 
mixture of particles with varying size ratio. Using a re- 
cursive approach, we obtained independent equations for 
the number of adsorbed particles of given sort and exact 
analytical expressions for the partial filling factors and 
variances for the larger particles. For the smaller parti- 
cles analytical expressions were obtained in a number of 
limiting cases. The results have been confirmed by direct 
Monte Carlo simulations. To do so, we have introduced 
a modified Monte Carlo procedure that enabled us to 
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explore a wide range of particle size ratios and particle 
fractions in the flux. 

A number of qualitative implications have been for- 
mulated, relevant to the energy branching problem in 
high-energy particle propagation through a semiconduc- 
tor crystal. Conclusions made concern the quantum yield 
and the energy resolution in semiconductor detectors 
made of crystals with several competing channels of im- 
pact ionization with different final electronic states. 

We have found a very strong anti-correlation effects 
which strongly suppress fluctuations of the total particle 
number compared to the fluctuations of partial contribu- 
tions by particles of a given sort. This effect is particu- 
larly evident when one considers the deposition of similar 
competing particles, e.g. parking of cars that arc differ- 
ent only in "color" . It may have dramatic consequences 
for semiconductor 7-radiation detectors, if the symmetry 
between anti-correlated particles is broken by a biased 
collection. This leads to an important conclusion that the 
energy resolution of semiconductor detectors is very sen- 
sitive to the collection efficiency of competing secondary 
particles. 

We have also found a very strong correlation effects 
that suppress fluctuations of the larger particle number 
for all particle ratios. As a result, the Fano factor for 
the larger particles is as a rule considerably smaller than 
that for the smaller particles. The variance of the cover- 
age by the smaller particles strongly increases with the 
growth of the particle size ratio b/a. This effect is due 
to the fluctuations in the size of gaps between larger par- 
ticles that serve as receptacles for small-particle deposi- 
tion. For b/a > 5 the small-particle variance exceeds that 
for the Poisson distribution in almost the entire range of 
particle fractions in the flux onto the adsorbing line. 
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APPENDIX A: SMALL PARTICLE 
CONTRIBUTIONS TO COVERAGE AND 
COVERAGE VARIANCE 



To calculate the contribution of small particles to the 
total coverage at large x, we use Eq. ^ with the initial 
boundary conditions With the substitution x — ^ x+b 
and using Eq. we rewrite Eq. ([T]) in the form 

{x + b — l)na{x + b) = q{b — a)na{b) + qx 

l*x+b—a px 

+2q na{y)dy + 2p na{y)dy (Al) 



Jb-a Jo 

Equation (jAip is valid for all x > b. Taking Laplace 
transformation of na{x) cut at a: < 6 by a step- function 
factor, we find that the transform, 

(A2) 



Na{s) = e ^''na{x)dx, 

Jb 



satisfies the following equation 



A + b-j] e'^Nais) = ±[l + {b - a)na{b)s] 
as I s- 



-t-2-e 
s 



Here 



'ljb-a)s 



Na{s)+Ji{s) 



2I 



Na{s) + J2{s) (A3) 



Ji{s) = I e~'"'na{x)dx, J2{s) ^ / e~'"'na{x)dx 



(A4) 



Rearranging the terms, we rewrite it in form 

1 



as s ^ ' 



(A5) 



where 



Rais)^q[l + ib-a)naib)s]+2s [qe^'''"'^' Jiis) + pJ2is) 

(A6) 
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The form of Eq. (|A5P is similar to Eq. ((20|) in which, 
however Ra should be calculated through Ji(s) and 
J2(.s). using Eqs. (|5l6p . For the case 6 < 2a we have 
na(i>) = 1, and Ji(s) ~ J2{s), while the explicit expres- 
sion for Ji(s) is easily obtained by substituting na(x) = 1 
in Eq. (|A4[) . Solution of Eq. (|A5|) then enables one to 
retrieve the result of Ref ^. To calculate Ji(s) and J2(s) 
for 6 > 2a, it is necessary to use Eq. ([5]), which describes 
RSA of small particles onto a short line a; < 6. Its ana- 
lytical solution and therefore the explicit expressions for 
Ji(s) and J2(s) can be obtained for the case b/a < 5 us- 
ing direct recursion to find na{x) (for one-particle RSA 
problem!). The result is rather cumbersome but suitable 
for numerical integration. 

For the case b/a > 5 one can exploit the exponentially 
rapid approach of the solution of Eq. ([6|) to its asymp- 



with 



totic behavior in the limit a; S> 1 (see e.g. 
numerical data). This asymptotic solution. 

R, 



271 for the 



riaix) 



(A7) 



can be used to calculate Ji(s) and then J2(s). To do this, 
we multiply Eq. ([6]) by exp(— sx) and integrate between 
and 6—1. We obtain an equation for J2(s) of the form 

J^{s) +(^a+ ^e— ^ Us) = -1 {e— /(s) 

+s{b - a)e-'"[naib) - 1] + 2se-"Vi(s)} , (A8) 



where 



r{b-a)s 

I{s) = I dyye~y. 



(A9) 

Solution of Eq. (jA8[) . satisfying the boundary conditions 
for Ua given by Eq. ([5]), is of the form 

•^2^ = 37TT^"" r dtm\naib)ib-a)te-^'-'^^' 
p{s)s'' Jo L 

+2tJi(t)- (l-e-"'-"'*)l (AlO) 



l3{t) = exp 



1 - 



dv 



(All) 



The contribution of small particles to the fill factor is 
then given by 



an. = I (i{u)Ra{u)du. 

10 



(A12) 



in which P{u) is given by the Eq. (|16[) and Ra{u) is 
defined by Eq. (|A6p . The obtained solution, though 
rather unwieldy, is suitable for numerical integration and 
for 6/a > 5 it gives the results that agree with Monte 
Carlo simulations. 

In the limiting case b/a ~ b/a ^ 1 it reduces to a more 
compact final expression for the contribution to the total 
coverage from the small particles 

6*^ = |l + ^ due-''"Pp{u) [q{u - 1) - 2pe"' 

(A13) 

with Pp{u) defined by For g = 1, Eq. (|XT3l) prop- 

erly gives 6a = R, while for q = one has 9a = R{1 — R). 
The latter expression corresponds to the coverage by 
small particles of the gaps between the large particles 
left after their initial deposition. For arbitrary q, the 
coverage given by Eq. (|A13p is depicted in Fig. 2 by the 
open squares. 

Similar approach can be used to calculate the small 
particle coverage variance. However, for b/a > 2 the 
equation for the Laplace transform of Ua{x) given by Eq. 
0. including all contributions to Na{s), becomes rather 
impractical. In the limiting case b/a ^ 1, when fiuctua- 
tions of the large particle gaps dominate the variance of 
small-particle coverage, one gets a more compact result 
shown in Fig. 5. 
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